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ABSTRACT 

Models of rapidly rotating massive stars at low metallicities show significantly 
different evolution and higher metal yields compared to non-rotating stars. 
We estimate the spin-down time-scale of rapidly rotating non-convective stars 
supporting an a — dynamo. The magnetic dynamo gives rise to mass loss in 
a magnetically controlled stellar wind and hence stellar spin down owing to 
loss of angular momentum. The dynamo is maintained by strong horizontal 
rotation-driven turbulence which dominates over the Parker instability. We 
calculate the spin-down time-scale and find that it could be relatively short, 
a small fraction of the main-sequence lifetime. The spin-down time-scale de- 
creases dramatically for higher surface rotations suggesting that rapid rotators 
may only exhibit such high surface velocities for a short time, only a small 
fraction of their main-sequence lifetime. 
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1 INTRODUCTION 

Models of rotating stars, particularly those of low-metallicity, have generated a great deal 
of interest because it is believed that their nucleosynthetic output could be significantly 
affected and so may better fit the observed composition of halo stars and globular clusters. 
Otherwise unaccountable nitrogen in very metal-poor stars could be explained by rapidly 
rotating massive stars because their nitrogen yields are found to greatly increase with spin 
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(.e.g. 



Meynet et al. 



20061 ) . There are suggestions that stars of very lo w metallicity cou 



very high equatorial spin velocities, of the order of 600 — 800 km s 1 ( IChiappini et al. 



d have 



2006J). 



Moreover very rapidly rotating stars could induce chemically homogeneous evolution and 
this may solve a proble m of insufficient mas s loss at very low metallicity in the collapsar 



mo dels of 7-ray bu r sts (jYoon &: Langerl 120051 ) . 



Tout &: Pringld ( 19911 ) showed how a magnetic dynamo might be driven by rapid rotation 



and radial turbulence. The mass loss in a magnetically controlled stellar wind results in 
quite rapid angular momentum loss and hence stellar spin down. In this paper we present 
an estimate of the spin-down rate for non-convective stars. In all cases the Alfven radius 
is larger than the stellar radius so angular momentum may be efficiently removed from the 
system and t his leads to short spin-down t i me-sc ales. 



Recently, 



Dervisoglu. Tout fc Ibanoglul (120101 ) considered the spin angular momentum 



evolution of the accreting components of Algol-type binary stars and demonstrated how 
accretion from a disc ought to spin the accreting star up to breakup velocity very quickly, 
after only a small amount of mass transfer. However, the accreting stars in Algols, which are 
hot main-sequence stars with radiative envelopes are not spinning at anywhere near their 
breakup rate. They showed that tides alone are insufficient to explain the spin down and 
argued that the most likely mechanism for the angular momentum loss is a magnetically 
controlled wind generated by a dynamo which is itself driven by the rapid rotation. This is 
further evidence that stars cannot continue to spin close to their breakup velocity and that 
a magnetic dynamo is a likely the spin-down mechanism. 



2 SPIN-DOWN TIME-SCALE AND MAGNETIC WIND 



We use as a basis the model of 



Tout fc Pringle! (11991 



hereinafter TP). In their model mass 



loss and braking are powered by an a — Q dynamo that extracts energy from differential 
rotation. The differential rotation at the surface and core boundary is maintained by stel- 
lar winds and core contraction respectively. Differential rotation throughout the radiative 
envelope results from redistribution of angular momentum because of the meridional cir- 
culation. This can be seen most clearly in models where angular momentum transport is 
treat ed as a purely diffusive proces s. In these models there is very little differential rota- 



tion (jPotter. Tout fc Eldridge 



201ll ). In t he case where there is a s trong magnetic field, the 



meridional circulation may be inhibited (IMaeder fc Meynetl 120031 ) but this conclusion has 
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been questioned by 



Maeder fc Meynet 



(120051 ) who found that, if the meridional circulation 
is neglected, there is insufficient differential rotation to maintain the magnetic dynamo. 
In particular, if differential rotation is too strongly suppressed by magnetic fields then it 
i s difficult to expla i n the existence of s low rotating stars with strong nitrogen enrichment 



(IHunter et al 



2009 



Meynet et al. 



2010l ). Mixing may be driven either by Parker's instability 



or by the horizontal rotation-driven turbulence that is itself responsible for the rotationally 
driven mixing. In TP's model convection served this purpose. Otherwise our a — Q dynamo 
operates in a similar way. Here we find that the extremely strong horizontal turbulence 
dominates over the Parker instability, even for low rotation rates. 

Our calculations are similar to those of TP. For a star of mass M, radius R* and angular 
velocity Q, the spin-down time-scale r S( j is the total angular momentum divided by the rate 
of angular momentum loss 
k 2 MR 2 Q 



T~sd 



(1) 



where k is the dimensionless stellar radius of gyration {k 2 ps 0.1) and J w is the angular 
momentum loss rate in the wind. 



J w « M w R 2 A tt, 



(2) 



where Ra is the Alfven radius ( IMestell Il968l ) where the magnetic energy and kinetic energy 
in the plasma balance. At Ra 



(3) 



where p w is the density of wind material and B is the magnetic flux density. We expect the 
velocity of the wind v w is of the order of the escape velocity so 



(2GM/R,) 



(4) 



We are imagining a wind launched by magnetic heating of a hot corona and accelerated 
by coronal pressure here. This differs from the radiatively driven winds normally associated 
with massive stars that would be somewhat faster. We have also assumed that the wind is 
not significantly slowed by the star's gravity. To quantify the effect of changing v w we write 

vl*P(2GM/R*), (5) 

where j5 > 1 if the wind is accelerated and < 1 if it is slowed. Though we expect Ra > R*, 
if it does fall below R* then the angular momentum loss rate is instead 



Jw 



M w Rtn 



(6) 
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because the material is then lost with the specific angular momentum of the stellar surface. 
Outside the star we assume the field strength B depends on radius as 

B = B*(RjR) n , (7) 

where B* is the field at the surface of the star. The outflow at the Alfven surface is roughly 
spherical so that 

p w « M w /(AjcR 2 a v w ) (8) 



and the relation between Ra and R* follows the equations (EJ) 



and dEb 



Ra 
R* 



( 



^BlRlV /{2n - 2) 



(9) 



\li M w v w J 

An undisturbed dipole field has n = 3. If the wind were able to open the field lines, so that 
n = 2, then Ra ^ R* would be larger and the spindown correspondingly faster. 

To estimate the wind mass-loss rate M w we consider the rate at which magnetic energy 
is dissipated and assume that all of this energy L w goes into launching the wind 

GMM W 



a- 



R* 



(10) 



We include the factor a because it is likely that not all the magnetic energy released is 
available to launch the wind. It might simply be radiated away. In that case a > 1. On the 
other hand it might even be possible that less energy is needed to drive the wind because it 
may be accelerated by radiation driving but only at high metallicity. In such a case a < 1. 

TP estimated the differential rotation by assuming that convection efficiently redis- 
tributes angular momentum. In non-accreting, non-convective, rotating stars the strength 
of the differential rotation results from the relative strength of meridional circulation and 
turbulent diffusion. The latter tends to smooth out differential rotation whereas the for- 
mer in general enhances it. Turbulent diffusion may arise owing to purely hydrodynamic or 
magnetic instabilities. The exact formulation with and without magnetic fields has been the 



subject of much discussion and many different models are common 



Heger et al. 



2000 



Meynet fc Maeder 



2000 



Maeder &: Meynet 



2004 



y used (for examples see 



Meynet et al 



2010). In 



all cases, the energy for driving the dynamo comes essentially from the excess rotational 
kinetic energy in the shear, -E s h ea r — ^k 2 MRlAQ 2 , where AQ approximates the range of 
angular velocity across the stellar radius. The rate L + at which energy is fed into the shear 
is 



1 k 2 MR 2 \9? 



(11) 
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R? i 

where t u fa -j- is the viscous time scale and v, which behaves as a viscosity, is the radial 
diffusion coefficient for angular momentum. 

We are considering radial variations of the rotation rate under a shellular framework 



(jZahnlll992l ) so it is appropriate to take the radial diffusio n coefficient for Kelvin-Hel mholtz 



Mevnet & Maederl (12000 ). Thus 



instabilities driven by the shear, referred to as -D s hear by 
AVl fa jrr^-\m is a measure of differential rotation in the star. Because -M^- varies throughout 
the star we choose r to give a typical value within the radiative envelope. We then evaluate 
the diffusion coefficients at the same radius. For massive stars, the amount of differential 
rotation is governed by the rate of evolution and the strength of the meridional circulation 



driven by rotation. Assuming that hydrostatic evolution is s 
and diffusion time-scales for angular momentum we may use 
state to get 

run 



ow compare d to the advection 



Zahnl (119921 ) 's model in steady 



AVI 



(12) 



where UP2(cos6) is the radial component of the merid ional circulation and P2{x) is the 



second Legendre polynomial in x. Using the formulation of 



Talon fc Zahnl ( 



visco sity which dominates the behaviour of the radial differential rotation (IMeynet fc Maeder 



1997) for the shear 



20001 ) we can approximate this as 
Att 



{2rUnN%\ 3 



\ K 



(13) 



where is the Brunt-Vaisala frequency and K is the thermal dif fusivity. It h as been 



show n that the stationary app roximation to the dynamical equation of I Zahnl (119921 ) is quite 
poor (IMeynet fc Maederl l2000l ). However this is in the case where the meridional circulation 
is inferred from the steady state differential rotation. Here we are using the meridional 
circulation and turbulent diffusivity calculated from non-stationary models to parameterize 
the differential rotation so we expect the approximation to be much better than a purely 
stationary model. We assume that the magnetic field strength does not become strong enough 
to cause dynamical feedback on either the rotation or the meridional circulation. 

Assuming that the energy dissipated from the shear is fed into driving a stellar wind, 
L + fa L w , we may derive the spin-down time-scales 



k 2 MR 2 n 



M w [max(R*,R A )} 2 Q 



_MoM 
AirBiR, 



i?A > R* 
Ra < R*- 



(14) 



The driver of the stellar dynamo is the shear, AQ, which converts the poloidal field B p to 
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toroidal Cyclon ic turbulence can then regenerate poloidal flux from toroidal (IParker 



1955 



Cowling! Il98lh . so we have 



~df 

and 

dB p 
~dT 



A Q Bp - ?± 

T 4> 



_r_ 



)B d 



B r 



(15) 



(16) 



where T is the dynamo regeneration term and r p and are the time-scales on which poloidal 
and toroidal flux are lost or destroyed. Following TP we take r p ~ « IOta, where 
ta ~ R*/va is the Alfven-wave crossing time in the star for an Alfven speed v\. As did TP, 
we write 



7*>t, 



(17) 



where 7 is an unknown parameter which describes the efficiency of the regeneration term 
and f t is a characteristic turbulent-ed dy velocity. Th e factor 7 is very uncertain. Based on 
the evolution of cataclysmic variables iwarner 1995 ) TP set it to 10 -2 for convective stars 



but foun d a much lower value of 3 x 10~ 5 to be sufficient to drive the activity in Ae and 
Be stars (ITout fc Pringldll994l ). We use a conservative estimate of 7 = 10~ 4 here and later 
show that our estimated spin-down rate is proportional to 7. 
In a steady state equations (|T5|) and (|T6l) lead to 

I'A 



B p « -^B, 



For each of the models calculated Ba 



0.1- 



-Ba. 



WB n so that 



Ba, W V A y/ floP*, 

where 

= M/(/i ^/ 3 ) 

is the mean density of the star. Equation (|T8|) then gives 
0.1 v\^pl 



B, 



and 



100.R,rAft. 



(18) 



(19) 



(20) 



(21) 



(22) 



We approximate the eddy velocity of turbulence owing to shear by v t ~ max(f !/ , O.lf a) 
where v v is the characteristic velocity of rotation-driven eddies. So either the rotation driven 
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turbulence dominates the dynamo regeneration or, if it is weak enough, the Parker instability 
sets in. TP set v*l ~ where is the horizontal diffusion coefficient. When rotationally 

driven turbulence dominates ^> v and is principally responsible for driving the dynamo. 
We find this to be the case in all models we consider below and so do not consider the 
alternative in any detail. 

It is the poloidal component of the magnetic field that predominantly drags around the 
wind so equations (TTT]) and ( 1201) - ( 1221) lead to 



B*^B P ^ 10^/3^0/4777^* W^ t . (23) 
We define time-scales 

r„ = ^, (24) 
v 

Tshear = AfT 1 (25) 

and 

^ = (Si) 2 (26) 



\2GM / 
so that 

M. = J^ (27) 

T shcar a T v 

and 

v w = P 1/2 ^- (28) 

7~dyn 

then with 

nurb = — • (29) 
Vt 

We then have the ratio of the Alfven to the stellar radius, 



Ra 



i 



™\l^i^liP^) , (30) 



Rt ~ \ V 47T fcl/2^V8 Uurb^dyn 

and a spin-down time-scale, 



gi/2gi/4fe r turb r shcar 7V R A> f> 



107-^/3/^0/4^ r 3/2 



dyn (31) 



r 2 r 

shear f D ^ D 

2 -TtA < Jt: 

r dyn 
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Table 1. Properties of a 20 Mq star spinning with equatorial velocity v (first column). The mass loss rate M/Movr - 1 is 
the predicted mass loss rate driven by the magnetic wind. The second column lists the ratio of the Alfven radius to the stellar 
radius. This ratio is proportional to a 1/4 /3— 1/8 7 1//2 . The star's surface magnetic field strength is B, and typical radial and 
horizontal diffusion coefficients are v and D^. Then AO is an estimate of the differential rotation through the star, r ms is its 
main-sequence lifetime, were it to continue to conserve its angular momentum, and r ac j is its spin-down time-scale which is 
proportional to a 1 ' 2 /3 1 ' 4 7 -1 . 



»/kms 1 


Ra/R* 


B*/G 


v 1 cm 2 s 1 


D h /cm 2 s- 1 


M/M Q yr 


-l 


ASl/s- 


i 


Tms/yr 


Wyr 


50 


2.1 


14 


1.0 x 10 8 


4.5 x 10 12 


2.7 x 10- 


10 


2.3 x 10 


— 5 


7.8 x 10 6 


1.7 x 10 9 


100 


1.7 


37 


3.3 x 10 s 


1.4 x 10 13 


3.9 x 10" 


9 


4.8 x 10 


-5 


7.9 x 10 e 


1.8 x 10 8 


150 


1.5 


62 


7.0 x 10 s 


2.5 x 10 13 


1.9 x 10" 


8 


7.4 x 10 


-5 


8.1 x 10 e 


4.6 x 10 7 


200 


1.4 


97 


1.0 x 10 9 


4.2 x 10 13 


5.3 x 10~ 


8 


9.9 x 10 


-5 


8.3 x 10 6 


1.8 x 10 7 


300 


1.3 


160 


2.0 x 10 9 


6.6 x 10 13 


2.2 x 10- 


7 


1.5 x 10 


-4 


8.9 x 10 6 


5.6 x 10 6 


400 


1.2 


220 


2.6 x 10 9 


9.4 x 10 13 


5.1 x 10- 


7 


2.0 x 10 


-4 


9.4 x 10 6 


2.7 x 10 6 


500 


1.2 


281 


3.0 x 10 9 


1.2 x 10 14 


9.0 x 10- 


7 


2.5 x 10 


-4 


9.7 x 10° 


1.6 x 10 6 



3 APPLICATION TO ROTATING MASSIVE STARS 



To investigate the spin down of a 20 M & star we evaluated the spin-down time-scale (equa- 
tion EI]), with a = 1, (3 = 1 and 7 = 10~ 4 , for a series of surface rotation rates. We used the 



Cambridge STARS code (lEggleton 
thors (most importantly by 



Pols et a] 



19711). which has previously been 



1995 



Stancliffe fc Eldridge 



updated by many au- 



2009). The version we have 



used has been modified to include the effects of rotation by lPotter. Tout fc Eldridgd (120111 ). 
With this we can estimate U, _D n and v for the star on the main sequence. In this case the 
Alfven radius is always larger than the radius of the star. We chose to estimate the turbulence 
at r where V — V a d has a minimum. In all cases r ~ 4 R Q from the centre of the star which 
has a total radius R* ~ 6R Q . In Tabled] we list various properties of the models, including 
the spin-down time-scale. We note that Ra/R* oc a 1 / 4 /3~ 1 / 8 7 1 / 2 and r S( j oc a 1 / 2 /? 1 ^" 1 . The 
mass-loss rate driven by this mechanism is at a rate of M w w a -1 10 -6 M Q yr" 1 for the 
most rapid rotators and escapes at v w « ^^kms" 1 . This is a slower, but not much 
slower, wind than that typically driven by radiation, at about 3,000 km s . So we may 
expect (5 ^ 3 and this to have little effect. Also there is no a priori reason why such a 
magnetically driven wind, or its consequences, should depend on metallicity. The spin-down 
time-scale becomes shorter than the main-sequence lifetime for stars initially rotating faster 
than around 200kms _1 . It is much shorter for rapid rotators and so a star cannot con- 
tinue to rotate at a high velocity throughout its entire main-sequence lifetime. It appears 
unlikely that homogeneous evolution, induced by very rapid rotation, can last long enough 
to change the fate of the stars even when they are of low metallicity. We note that ( Pup- 



pis, a 53.9M star flRepolust. Puis fc Herrero 
v sini 200kms _1 ( Kudritzki. Simon fcHamann 



200 41) wi th an observed rotational velocity of 



19831 ). has an observed mass loss rate of 
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about M w m 3.5 x 10 6 M yr 1 (j Cohen et al.ll2010l ). This mass loss rate is not inconsistent 
with our estimate. 



3.1 Uncertainty in the estimation 

We have assumed perfect efficiency for the transfer of energy from shear to the magnetically 
driven wind. This is unlikely to be the case and so the spin-down time-scales are likely to 
be larger than those predicted. However, if the efficiency factor a does not depend on the 
rotation rate, the prediction that rapid rotators spin down on a much shorter time-scale is 
still valid. Indeed, unless the efficiency is extremely low (a ^ 30), the spin-down time-scale 
for rapid rotators is much shorter than their main-sequence lifetime. 

There is also uncertainty introduced by our choice of vq. We have chosen it by a stan- 
dardised method within the radiative envelope. However, if we vary ro for a 20 M star 
with surface rotation of 200 km s" 1 , the spin-down time-scale varies between 7.2 x 10 4 and 
2.1 x 10 7 yr. It is unlikely to be close to the lower limit because this is the spin-down time- 
scale when r is taken close to the convective core. This is only representative of a small 
part of the radiative envelope and energy deposited there is unlikely to be very influential 
in driving the stellar wind. 

A major uncertainty is 7, the standard dynamo regeneration term. In our calculation, 
we have used a conservative value of 10 -4 but it could be larger by as much as two orders of 
magnitude. If that were the case then, the spin-down time-scales (r s d oc 7 _1 ) would be much 
shorter. Given the uncertainties in our estimate, we cannot draw a quantitative conclusion 
for a star with a moderate spin rate. However, our results do suggest that stars rotating near 
their break-up velocity, or stars with very efficient rotationally driven mixing, are likely to 
spin down very quickly by driving a magnetically locked wind. This effect cannot be ignored. 

As table [1] shows, Ra > R* in all of the cases examined here. This means that r S( j oc AQ, 
assuming that the turbulent viscosity is unaffected. If the differential rotation is strongly 
suppressed by the presence of a magnetic field then it may extend the spindown timescale 
of the star. However, the degree to which this happens is extremely uncertain. In partic- 
ular, if solid body r otation is enforced it is difficult to reconcile models and observations 
(IMevnet et al.lboioh . 
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4 CONCLUSION 

We have estimated the spin-down time-scale for rapidly rotating non-convective stars. With 
a conservative value for the magnetic field regeneration term 7, the spin-down time-scale 
is likely to be long for slowly rotating stars for which the horizontal turbulence driven by 
rotation is weaker. However, for more rapid rotators, the spin-down time-scale can be much 
shorter than the main-sequence lifetime. Given the uncertainties in the coefficient of mixing, 
viscosity and the appropriate mixing length, more detailed models of magnetic spin down 
driven by rotational mixing must be developed before we can decide on whether rapid spin 
can lead to homogeneous mixing and on the importance of rotationally driven mixing at 
lower spin rates. These calculations do however mean that we should not automatically 
assume that even very low-metallicity stars can continue to spin rapidly throughout their 
main- sequence lives. 
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